Casimir interactions of an object inside a spherical metal shell 
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We investigate the electromagnetic Casimir interactions of an object contained within an otherwise 
empty, perfectly conducting spherical shell. For a small object we present analytical calculations of 
the force, which is directed away from the center of the cavity, and the torque, which tends to align 
the object opposite to the preferred alignment outside the cavity. For a perfectly conducting sphere 
as the interior object, we compute the corrections to the proximity force approximation (PFA) 
numerically. In both cases the results for the interior configuration match smoothly onto those for 
the corresponding exterior configuration. 



The development of methods to measure Casimir forces 
with high precision |lj and the advent of microelectrome- 
chanical devices, for which Casimir forces are impor- 
tant |21 E] , have stimulated the development of new the- 
oretical tools for analyzing these forces for geometries 
beyond parallel plates [4j. In particular, methods based 
on scattering theory O |6l 1^ have been used to precisely 
compute Casimir forces and torques for a wide variety of 
configurations of both perfect conductors and dielectrics. 
We have applied these methods to find the electrody- 
namic Casimir interaction of a conducting or dielectric 
object inside a perfectly conducting spherical cavity [Sj. 
Earlier studies of Casimir interactions of one object con- 
tained within another have been limited to infinitely long 
cylinders 

In this paper we report two sets of results. First, 
we determine the forces and torques on a small object, 
dielectric or conducting, well separated from the cav- 
ity walls. This is the interior analogue of the famous 
Casimir-Polder force on a polarizable molecule near a 
perfectly conducting plate [TT]. Second, we compute the 
interaction energy of a finite-size metal sphere with the 
cavity walls when the separation, d, between their sur- 
faces tends to zero. The leading term in d is known to be 
given by the PFA . By evaluating our closed- form ex- 
pressions numerically, we find the next to leading term in 
an expansion in d. The main achievements of our analy- 
sis are: (i) The first exact results for the electromagnetic 
interaction of an interior macroscopic object with a con- 
fining cavity; (ii) The determination of the corrections 
to the PFA for two spheres for all values of their radii, 
both positive (exterior) and negative (interior); (iii) A 
physical understanding of the absence of an orientation 
dependence of the Casimir-Polder (CP) force between an 



object and a plane mirror; and (iv) The generalization of 
the CP force to concave mirrors. The last follows from 
our analysis when the cavity radius is much larger than 
d and the size of the interior object. Then the interac- 
tion is determined only by a concave segment of the shell 
that is closest to the object and our results are expected 
to apply also to this open geometry. Our work suggests 
that the orientation dependence of Casimir forces near 
curved surfaces might prove to be an interesting area for 
future experiments. Furthermore, the interactions de- 
scribed here are relevant for trapped atoms in spherical 
optical resonators ^13j . 

An understanding of the corrections to the PFA has 
been sought for some time T3]. Our works extends be- 
yond previous attempts, many of which have treated hy- 
pothetical scalar fields, by considering electromagnetic 
fields for two objects with different curvature and rela- 
tive position. For the sphere-plane geometry jjlSlITH] and 
two spheres of equal radii facing each other [5] the cor- 
rections have recently been computed. We repeat these 
computations for other ratios of radii of two spheres out- 
side one another and find that the inside and outside 
results connect smoothly. This is of direct experimen- 
tal significance because Casimir force measurements are 
usually performed with spherical rather than perfectly 
planar surfaces. 

Our analysis is facilitated by a recently developed for- 
malism that combines path integral and scattering theory 
techniques and is applicable to general shapes and config- 
urations, including objects inside each other [HIISIIIZ]- In 
the case where an object, S^, lies inside a perfectly con- 
ducting cavity with inner surface Eq, the Casimir energy 
is given by [3 [17] , 
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J^i is the scattering amplitude for electromagnetic waves 
off the interior object; JFj"^ is the inverse scattering am- 
plitude for the conducting cavity, a sphere in our case. 
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The appearance of the inverse of To is a result of the 
interior geometry |17j . These scattering amplitudes are 
matrices evaluated in a spherical vector wave basis with 
respect to appropriately chosen origins within each ob- 
ject. The translation matrices, and V*°, relate regu- 
lar wave functions between the coordinate systems of the 
interior object and the spherical cavity. (yV"° ~ V*°^ up 
to multiplication by (—1) of some matrix elements, see 
Ref. [T7] for details.) All of these matrices are functions 
of the imaginary frequency uj = ik. The determinant in 
the denominator subtracts the Casimir energy when the 
origins of the two objects coincide. This way of normal- 
izing the Casimir energy differs from the exterior case 
where the objects are removed to infinite separation; a 
choice that would be unnatural in the interior case. 

Let us first consider the limit, in which the interior 
object is much smaller than the radius of the spher- 
ical cavity. Using the matrix identity IndetA^ = 
TrlnA^, the integrand in Eq. (fTl) can be expanded as, 
E = -hc/{2TT) dKTv{M + + ...), where N = 

TiV^" describes a reflection event where a wave 
travels from one object to the other and back 5 . In gen- 
eral, all terms in the series expansion are important, il- 
lustrating the fundamentally non-two-body nature of the 
Casimir force. The rate of convergence of this series de- 
pends on the size of relative to the separation of its 
surface from that of Sq. 

If the inner object is small compared to the size of 
the cavity, the first term in the series, Tr A/", already 
gives an excellent approximation to the energy. In the 
small-size limit, the scattering amplitude Ti^imP,i'm'P' i 
(where I and m are angular momentum indices and P la- 
bels M or E polarization) can be expanded in powers of 
K. Only the following terms contribute to lowest order: 
Ti^unP.Un'pin) = 2K^a^^,/3 + 0{K'^), where a^„, is the 
static electric (P = E) or magnetic (P — M) polarizabil- 
ity of the inner object, the same tensor that determines 
the Casimir- Polder interaction [TTJ [T51 [12] . 

We fix So to be a spherical shell of radius R and define 
a to be the displacement of the center of the interior 
object from the center of the sphere. We find for the 
Casimir energy to leading order in r/R (where r is the 
typical length scale of the interior object), the interior 
analog of the Casimir-Poldcr interaction ^ , 
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The z-axis is oriented from the center of So to S^ , and 
represent the interior object's static polarizability tensors 
in a Cartesian basis. Corrections to this energy come 
from dynamic (frequency dependent) dipole polarizabil- 
ities and higher order multipoles |20l . The coefficient 
functions and , plotted in Fig.nJ can be expressed 
in terms of modified Bessel functions and Ki, as (with 
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FIG. 1: Plot of the functions f^^^'^{a/R) and g"/'^(a/i?), 
defined in Eqs. (|3| and Q, respectively. 
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and / and g are obtained by substituting (P ^ M) 



in the above equations. The functions 

(with ■q = l + \/2) 



M/E 



are given by 



_ K,{x) . _ K,{x)+2xK'^{x) 
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is negative and decreasing with a/R, while /^^ is 
positive and increasing. Notice that the g^ are about 
an order of magnitude smaller than the . Therefore, 
it is safe to conclude that an interior object will always 
be attracted to the spherical cavity walls if <C . 
If one could make a material, for which a*^ was large 
and dominant, the interior object could float inside the 
spherical shell instead of being attracted to its walls. 

There are important differences between Eq. ^ and 
the classic Casimir-Polder result: first, the energy de- 
pends in a non-trivial way on a/R; second, at any non- 
zero distance from the center, the interior object expe- 
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riences a torque; and third, the force between the two 
bodies depends on the interior object's orientation. 

To explore the orientation dependence of Eq. (|2| as- 
sume, for simplicity, there is a single frame in which both 
and a^^ are diagonal. In this body-fixed frame, write 
- «TO = /3 and - ^{a%^ + a^y) = j (where we 
have suppressed the M/E label). The polarizability in 
the "lab frame" is obtained by a = TZa^TZ"^ , where TZ 
is a rotation matrix that orients the principal axes of the 
inner object with respect to the lab frame [19]. This pro- 
cedure leaves Tra*^ invariant, and gives for the second 
Hne in Eq. ([2|, 



9^{a/R) ( sin^ 6 cos 20 + 7-f'(3 cos^ 9-1) 

P=M,E 



where 4> corresponds to the azimuthal rotation of the ob- 
ject about its principal z-axis, and 6 is the angle between 
the object's principal z-axis and the "laboratory" z-axis 
connecting the center of the sphere to the origin of S^. 

If /3 then the object held at fixed inclination, 6*, 
experiences a torque that causes it to rotate about the 
body-fixed z-axis. If, however, the object has axial sym- 
metry (/3 = 0), then the only torque on the object tries to 
align it either parallel or perpendicular to the displace- 
ment axis. 

A "cigar shaped" object (7 > 0) prefers to orient so 
as to point perpendicular to the z axis, and a "pancake" 
(7 < 0) tries to align its two large axes perpendicular to 
the z axis. The small ellipse inside the sphere in Fig.[2]il- 
lustrates a side view of both the cigar and the pancake in 
their preferred orientation. It is interesting to note that 
and are both positive. So, in contrast to the force, 
the contributions to the torque from magnetic and elec- 
tric polarizabilities are in the same direction, if they have 
the same sign. More complicated behavior is possible if, 
for example, the electric and magnetic polarizabilities are 
not diagonal in the same body-fixed coordinate system. 
Note that our results cannot be compared to the PFA 
approximation since the the size of the inner object, not 
the separation of surfaces, d, has been assumed to be the 
smallest scale in the analysis. 

An identical analysis can be performed for a polariz- 
able object outside a metallic sphere where a/R > 1. The 
analogous exterior functions f{a/R) and g{a/R) are ob- 
tained by exchanging I^, and K^, in Eqs. Q and 
It turns out that g{a/R) < for both polarizations. 
Therefore, the preferred orientation of a polarizable ob- 
ject outside a metallic sphere is opposite of that in the 
interior case (see the small ellipse outside the large sphere 
in Fig. [2]). The continuation of the functions / and g from 
"interior" to "exterior" is displayed in Fig. [2] where the 
transition from one orientation to the other is clear. 

In the limit x = a/R ^ 1, we find that the func- 
tions /^/^(x) diverge like (1 - a;)-^ while g^^'^'ix) di- 
verge only as {l — x)~^. These reproduce the well-known 
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FIG. 2: The ratio 3^//^, which dete rmines the preferred 
orientation of the interior object, plotted versus x — a/R 
showing the change in preferred orientation from interior 
{a/R < 1) to exterior {a/R > 1) (displayed by two small 
ellipses as described in the text). The data points are numer- 
ically computed from Eqs. ([sjl, H. The solid curves are fits 
of the form ci(l — x) + C2{1 — x) to these data points. 



Casimir-Polder results for an object facing a plane: the 
object's energy grows like 1/d* and does not depend on 
its orientation at leading order. 

The results, presented up to this point, are accurate 
in the limit r/R — > 0, where r is the typical length scale 
of the interior object. For fixed r/R they become more 
accurate as the separation from the center, a/R, vanishes. 
For example, for a metallic sphere of radius r, at a/R = 
0.2, Eq. ([2| differs from the exact result by 5% for r/R = 
0.1 and 20% for r/R = 0.2. A more detailed comparison 
can be found in Ref . [8, . 

The opposite extreme from the small-object limit 
comes when the interior object is nearly touching the 
cavity wall. In this limit the Casimir force F between 
two conducting spheres, which is attractive, is propor- 
tional in magnitude to d~^, where d — R — r — ais the 
separation of surfaces. The coefficient of is given by 
the PFA [12], 



Km d^F^ 



TT^hc rR 
360 r + R 



(6) 



This result holds for both the interior and the exterior 
configuration of two spheres. For fixed r we formally 
distinguish the cases: R > for the exterior, i? — > cx) for 
the plate-sphere, and R < for the interior configuration 
(see Fig. [3] for reference). All possible configurations are 
taken into account by considering — l<r/i?<l. 

Recently, there has been much interest in determin- 
ing the first correction to the leading PFA result to order 
d/r [131 [TS]. The plate-sphere case and exterior problems 
of spheres of equal radii have been computed in Refs. [S] 
and ^15j , respectively. Since most experiments up to now 
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have considered spherical conductors separated by dis- 
tances much smaller than their radii, the first correction 
in d/r to the PFA is the geometric correction of greatest 
practical interest. 

Our numerical studies of the interior configuration 
for conducting spheres [5] enable us to study the limit 
d/r — > 0. This is a difficult limit because no simplifying 
approximations can be applied to Eq. ([ij. All powers 
of J\f contribute, and the number of partial waves (/, V) 
necessary to obtain convergence grows as d/r — s- 0. Al- 
though we know of no derivation of the functional form 
of the Casimir force beyond the leading term in the PFA, 
our numerical data are well fit by a power series in d/r, 



TT^hc rR 
' 360d3 r + R 
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We have used this functional form to extract the coeffi- 
cient 9i{r/R). 

It is useful to have an estimate, however crude, of 
9i(r/R) over the whole range of r/R with which to com- 
pare our results. Although the PFA is accurate only in 
the limit d/r 0, it can be extended in various ways 
to the whole range of d, r, and R. The PFA is obtained 
by considering both surfaces as made up of infinitesimal 
parallel mirrors. From each point (^1,^2) on the surface 
of object O one computes the distance ^(^1,1^2) to the 
other object's surface along the surface normal n(^i,^2)- 
By integrating the Casimir energy per unit area for two 
parallel plates separated by L(^i,<^2) over the surface of 
object O one obtains the "O-based" PFA energy. Clearly, 
the result depends on which object one chooses as O, but 
the various results do agree to leading order in d/r. We 
can choose either of the two spheres to arrive at the "r- 
based PFA" or the "i?-based PFA" , see Fig. [S] Either 
one yields a 'correction' to the leading order PFA, 
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where x — r/R. Again, 9f^^ and are only used for 
comparison with the actual correction 9i. Note that the 
PFA predicts a smooth continuation from the interior to 
the exterior problem. 

In Fig. [3] we plot the values of 9i extracted from a 
numerical evaluation of the force from Eq. ([T]) for various 
values of r/R < 0, along with the values for r/R = and 
r/R = 1 from Refs. [5] and [15]. We have also repeated 
the exterior analysis of Ref. [5] for other values of r/R > 
0. For reference, the two PFA estimates are also shown. 

Eq. ([1]) is numerically evaluated by truncating the ma- 
trix M at finite multipole order I, and extrapolating to 
obtain the I 00 limit. For the data in Fig. |3] Af was 
truncated at / = < 60 for the interior and at Z = /' < 35 
for the exterior configuration. 

The numerical data in Fig. [3]show a smooth transition 
from the interior to the exterior configuration. Although 
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FIG. 3: PFA correction coefficients for spheres. r/R ranges 
from -1 (interior concentric), to zero (sphere-plane), to -1-1 
(exterior, equal radii). The data points correspond to the ex- 
act values of 61 calculated numerically, while the solid black 
curve is a fit (see text). Inset: "interior" and "exterior" geo- 
metrical configurations. 



the PFA estimates do not describe the data, the r-based 
PFA has a similar functional form and divergence as x — > 
— 1. Therefore, wc fit the data in Fig. |3]to a function, 
9i{x) = -{kix + k2x/{l + x) + ks) and find, fci = 1.05 ± 
0.14, k2 = 1.08 ± 0.08, fca = 1.38 ± 0.06. This provides a 
simple form for the leading PFA correction for metallic 
spheres, one inside the other and both outside, which is 
relevant for many experiments. Notice, however, that 
the actual function 9i{x) is not known analytically and 
that our fit represents a reasonable choice which may not 
be unique. Our results show that the correction to the 
PFA has a significant dependence on ratio of curvatures 
of the two surfaces. The correction is a factor of two 
larger for two spheres of equal radii than for the sphere- 
plane setup; it vanishes near r/R — —0.5; and it becomes 
large positive as r/R —1. These effects should be 
taken into account in future experimental searches for 
PFA corrections. 
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